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Abstract 

We discuss the role of absorbtion and antishadowing in particle produc- 
tion. We reproduce power-like energy behavior of the mean multiplicity in 
the model with antishadowing and discuss physical implications of such be- 
havior for the hadron structure. 
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Introduction 



Multiparticle production and global observables such as mean multiplicity and 
its energy dependence alongside with total, elastic and inelastic cross-sections 
provide us a clue to the mechanisms of confinement and hadronization. General 
principles are very important in the nonperturbative sector of QCD and unitarity 
which regulates the relative strength of elastic and inelastic processes is the one of 
such principles. There is no universal, generally accepted method to implement 
unitarity in high energy scattering and as a result of this fact a related problem 
of the absorptive corrections role and their sign has a long history (cf. [JT)] and 
references therein). However, a choice of particular unitarization scheme is not 
just a matter of taste. Long time ago arguments based on analytical properties 
of the scattering amplitude were put forward [Q| in favor of the rational form 
of unitarization. It was shown that this form of unitarization reproduced correct 
analytical properties of the scattering amplitude in the complex energy plane much 
easier compare to the exponential form. 

Interest in unitarity limitations and the respective dependencies of the global 
observables was stimulated by the preparation of the experimental program for 
the LHC and the future plans to study soft interactions at the highest available 
energies. For example, correct account for unitarity is essential under theoretical 
estimates of the Higgs production cross-section via diffractive mechanisms. The 
region of the LHC energies is the one where new, so called antishadow scattering 
mode can be observed. Such a mode naturally appears when the rational form of 
unitarization being exploited It has been demonstrated that this mode can be 
revealed at the LHC directly measuring a e i(s) and a to t(s) 0] (and not only by 
means of the analysis of impact parameter distributions). Antishadowing leads 
to self-damping of the inelastic channels and asymptotically dominating role of 
elastic scattering, i. e. cr e i(s) / a tot (s) — > 1 at s — > oo. Immediate question 
arises on consistency of this mechanism with the growth of mean multiplicity in 
hadronic collisions with energy. Moreover, many models and experimental data 
suggest a power-ilke energy dependence of mean multiplicity!] and a priori the 
compatibility of such dependence with antishadowing is not evident. 

In this note we apply a rational (U -matrix) unitarization method to con- 
sider the global features of multiparticle dynamics such as mean multiplicity and 
role of absorptive corrections. We show that it is possible to reproduce power-like 
energy behavior of the mean multiplicity in the model with antishadowing and 
discuss some physical implications. 

'Recent discussions of power-like energy dependence of the mean hadronic multiplicity and 
list of references to the older papers can be found in [pi 
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1 Multiparticle production in the U -matrix approach 



The rational form of unitarization is based on the relativistic generalization of the 
Heitler equation of radiation dumping [§]. In this approach the elastic scattering 
amplitude satisfies unitarity since it is a solution of the following equation 

F = U + iUDF (1) 

presented here in the operator form. Eq.[l] allows one to fulfill the unitarity pro- 
vided the inequality 

hnU(s, b)>0 (2) 
is satisfied. The form of the amplitude in the impact parameter representation is 

/M) = T i %^T Y , (3) 
1 — iU{s, b) 

where U (s, b) is the generalized reaction matrix. It is considered as an input dy- 
namical quantity similar to the eikonal function. It is to be noted that the analogous 
form for the scattering amplitude was obtained by Feynman in his parton model 
for diffractive scattering which he has never published (cf. [0]). 

In the impact parameter representation the unitarity equation rewritten at high 
energies for the elastic amplitude f(s, b) has the form 

lmf(s,b) = \f(s,b)\ 2 + V (s,b) (4) 

where the inelastic overlap function 

1 d<7i ne l 



rj(s,b) = 



4.7T db 2 



is the sum of all inelastic channel contributions. It can be expressed as a sum of 
n-particle production cross-sections at the given impact parameter 

V(s,b) = ^2a n (s,b), (5) 

n 

where 

1 d<j f°° 
a n (s,b) = — -jTT, crn(s) = 8tt / bdba n (s,b). 
4.7T db 2 J 

Inelastic overlap function is related to U(s,b) according to Eqs. (TJ) and (Q) as 
follows 



3 



Then the unitarity Eq. ^ points out that the elastic scattering amplitude at 
given impact parameter value is determined by the inelastic processes when the 
amplitude is a pure imaginary one. Eq. |] imply the constraint \f(s,b)\ < 1 
while the "black disk" limit presumes inequality \f(s, b)\ < 1/2 and the elastic 
amplitude satisfying the latter condition is a shadow of inelastic processes. In its 
turn the imaginary part of the generalized reaction matrix is the sum of inelastic 
channel contributions: 

ImU(s,b) = J2u n (s,b), (7) 

n 

where n runs over all the inelastic states and 

U n (s,b) = J dT n \U n (s,b,{^ n })\ 2 . (8) 

In Eq. (|8|) dT n is the n-particle element of the phase space volume. The functions 
U n (s, b, {£ n }) are determined by dynamics of hi + h 2 — > X n processes, where 
{£„} stands for the full set of respective kinematical variables. Thus, the quantity 
ImU (s, b) itself is a shadow of the inelastic processes. However, unitarity leads to 
self-damping of the inelastic channels pSp in the sense that increase of the func- 
tion Imf/(s, b) results in decrease of the inelastic overlap function r)(s, b) when 
Im[/(s, b) exceeds unity (cf. Fig. 1). 




b = R (s) b b=R(s) b 



Figure 1: Shadow and antishadow scattering regions 

Respective inclusive cross-section [KJ which takes into account unitarity 
in the direct channel has the form 

d£ J |1 - iU{s,b)\ 2 

The function I(s, b, £) in Eq. (g) is expressed via the functions U n (s ) 6, £, {Cn-i}) 
determined by the dynamics of the processes hi + h 2 — ► h 3 + X n _i. 



H s ' b iO = ^2 n / dr n \u n (s,b,£, {^ n _i 

n>3 



})| 2 (10) 
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and 

I(s,b,£)d£ = n(s,b)lmU(s,b). (11) 



The kinematical variables £ (x and p±, for example) refer to the produced 
particle h 3 and the set of variables {£ n -i} describe the system X n _i of n — 1 
particles. 

Now we turn to the mean multiplicity and consider first the corresponding 
quantity in the impact parameter representation. As it follows from the above the 
n-particle production cross-section a n (s, b) 



Then the probability 



a ^ b)= \l- i U(s,bW (12) 



p t u\- a n(s,b) U n (s,b) 

Pn{s, b) = — - = - — — — . (13) 

Vinei{s,b) JmU(s,b) 



Thus, we observe the cancellation of unitarity corrections in the ratio of the 
cross-sections a n (s, b) and cr ine ;(s, b). Therefore the mean multiplicity in the im- 
pact parameter representation 



b) = J2 nP n(s,b) 



is not affected by unitarity corrections and therefore cannot be proportional to 
rj(s, b). This conclusion is consistent with Eq. (|TT|). The above mentioned pro- 
portionality is a rather natural assumption in the framework of the geometrical 
models, but it is in conflict with the unitarization. Because of that the results 
[ |TT| ] based on such assumption should be taken with precaution. However, the 
above cancellation of unitarity corrections does not take place for the quantity 
n(s) which we address in the next section. 



2 Growth of mean multiplicity 

We use a model for the hadron scattering described in [|12[j. It is based on the 
ideas of chiral quark models. The picture of a hadron consisting of constituent 
quarks embedded into quark condensate implies that overlapping and interaction 
of peripheral clouds occur at the first stage of hadron interaction (Fig. 2). Non- 
linear field couplings could transform then the kinetic energy to internal energy 
and mechanism of such transformations was discussed by Heisenberg [|T3p and 
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Carruthers [|T^]. As a result massive virtual quarks appear in the overlapping re- 
gion and some effective field is generated. Valence constituent quarks located 
in the central part of hadrons are supposed to scatter simultaneously in a quasi- 
independent way by this effective field. 




Figure 2: Schematic view of initial stage of the hadron interaction. 

Massive virtual quarks play a role of scatterers for the valence quarks in elastic 
scattering and their hadronization leads to production of secondary particles in the 
central region. To estimate number of such scatterers one could assume that part 
of hadron energy carried by the outer condensate clouds is being released in the 
overlap region to generate massive quarks. Then this number can be estimated by: 

N( s b) oc (1 - ® = N (s)D c (b), (14) 

m Q 

where txiq - constituent quark mass, (ko) - average fraction of hadron energy car- 
ried by the constituent valence quarks. Function D\ describes condensate distri- 
bution inside the hadron h, and b is an impact parameter of the colliding hadrons. 

Thus, N(s, b) quarks appear in addition to N — n hl + n h2 valence quarks. 
In elastic scattering those quarks are transient ones: they are transformed back 
into the condensates of the final hadrons. Calculation of elastic scattering am- 
plitude has been performed in [JT2|] . However, valence quarks can excite a part 
of the cloud of the virtual massive quarks and these virtual massive quarks will 
subsequently fragment into the multiparticle final states. Such mechanism is re- 
sponsible for the particle production in the fragmentation region and should lead 
to strong correlations between secondary particles. It means that correlations exist 
between particles from the same (short-range correlations) and different clusters 
(long-range correlations) and, in particular, the forward-backward multiplicity 
correlations should be observed. This mechanism can be called as a correlated 
cluster production mechanism. Evidently, similar mechanism should be signif- 
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icantly reduced in e + e~-annihilation processes and therefore large correlations 
are not to be expected there. 

As it was already mentioned simple (not induced by interactions with valence 
quarks) hadronization of massive N(s, b) quarks leads to formation of the multi- 
particle final states, i.e. production of the secondary particles in the central region. 
The latter should not provide any correlations in the multiplicity distribution. 

Remarkably, existence of the massive quark-antiquark matter in the stage pre- 
ceding hadronization seems to be supported by the experimental data obtained at 



CERN SPS and RHIC (see [fj] and references therein). 

Since the quarks are constituent, it is natural to expect direct proportionality 
between a secondary particles multiplicity and number of virtual massive quarks 
appeared (due to both mechanisms of multiparticle production) in collision of the 
hadrons with given impact parameter: 

n{s, b) = a(n hl + n h2 )N (s)D F (b) + f3N (s)D c (b), (15) 

with constant factors a and (3 and 

D F (b) = D Q ® D c , 

where the function Dq(b) is the probability amplitude of the interaction of valence 
quark with the excitation of the effective field, which is in fact related to the quark 
matter distribution in this hadron-like object called the valence constituent quark 



QI2]. The mean multiplicity n(s) can be calculated according to the formula 



= 1Z ' n ( s Mv{s,b)bdb (i6) 
rj(s,b)bdb 

It is evident from Eq. (|16|) and Fig. 1 that the antishadow mode with the periph- 
eral profile of i](s, b) suppresses the region of small impact parameters the main 
contribution to the mean multiplicity is due to peripheral region of b ~ R(s). 

To make explicit calculations we model for simplicity the condensate distri- 
bution D c (b) and the impact parameter dependence of the probability amplitude 
Dq(b) of the interaction of valence quark with the excitation of the effective field 
by the exponential forms, and thus we use exponential dependencies for the func- 
tions Dp{b) and Dc(b) with the different radii. Then the mean multiplicity 

n{s, b) = aN (s) exp(-b/R F ) + (3N (s) exp(-b/R c ). (17) 

The function U (s, b) is chosen as a product of the averaged quark amplitudes 

N 

U(s,b) = H(f Q (s,b)) (18) 

Q=l 
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Figure 3: Energy dependence of mean multiplicity, theoretical curve is given by 
the equation n(s) = as s (a = 2.328, 5 = 0.201); experimental data from the Refs. 

[in. 



in accordance since in the model valence quarks scatter in effective field simulta- 
neously and quasi-independently. The 6-dependence of (Jq) related to the quark 
formfactor Fg(q) has a simple form (fq) oc exp(— mg&/£). Thus, the generalized 
reaction matrix (in a pure imaginary case) gets the following form [|T2"|] 



U(s,b) = ig 



1 + a 



m Q 



N 



eM-Mb/0, 



(19) 



where M = ^2^ =1 mQ. At sufficiently high energies where increase of the total 
cross-sections is quite prominent we can neglect the energy independent term and 
rewrite the expression for U(s, b) as 



U{a,b) = ig {s/m Q ) N/ \M-Mb/i). 



(20) 



After calculation of the integrals ( [R)| ) we arrive to the power-like dependence 
of the mean multiplicity n(s) at high energies 



n(s) 



as 



s F 



bs 5c , 



(21) 



where 



; i- 



ITLqR! 



and 5 f 



c 



; l- 



m Q R 



c 



There are four free parameters in the model, a, (3 and Rf, Rc, and the freedom 
in their choice is translated to a, b and 8p, 5c- The value of £ = 2 is fixed from 
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the data on angular distributions [ Jl"2"l ] and for the mass of constituent quark the 
standard value mg = 0.35 GeV was taken. However, fit to experimental data 
on the mean multiplicity leads to approximate equality 5p — 5c and actually 
Eq. ( PH ) is reduced to the two-parametric power-like energy dependence of mean 
multiplicity 

n = as s , 

which is in good agreement with the experimental data (Fig. 3). Equality 5p — 5c 
means that variation of the correlation strength with energy is weaker than the 
power dependence and could be, e.g. a logarithmic one. From the comparison 
with the data on mean multiplicity we obtain that 5 ~ 0.2, which corresponds to 
the effective masses, which are determined by the respective radii (M = l/R), 
M c ~ M F ~ 0.3m Q , i.e. Mp ~ Mc — m n . 

The value of mean multiplicity expected at the LHC energy (y/s = 14 TeV) 
is about 1 10. It is not surprising that it is impossible to differentiate contributions 
from the two mechanisms of particle production at the level of mean multiplicity. 
The studies of correlations are necessary for that purpose. 

Conclusion 

It was shown that the model based on accounting unitarity [ |T2"| ] and extended to 
multiparticle production provides a reasonable description of the energy depen- 
dence of mean multiplicity leading to its power-like growth with a small expo- 
nent. This result is a combined effect of unitarity and existence of the phase 
preceding hadronization when massive quark-antiquark pairs are generated. It is 
worth noting again that power-like energy dependence of mean multiplicity ap- 
pears in various models and is in a good agreement with heavy-ion experimental 
data tooQ. 

Multiplicity distribution P n (s, b) and mean multiplicity n(s, b) in the impact 
parameter representation have no absorptive corrections, but since antishadow- 
ing leads to suppression of particle production at small impact parameters and 
the main contribution to the integral multiplicity n(s) comes from the region of 
b ~ R(s). Of course, this prediction is to be valid for the energy range where 
antishadow scattering mode starts to develop (the quantitative analysis of the ex- 
perimental data [[18]] gives the value: ^/so ~ 2 TeV) and is therefore consistent 
with the "centrality" dependence of the mean multiplicity observed at RHIC [Jl~9|] . 

It is also worth noting that no limitations follow from the general principles for 
the mean multiplicity, besides the well known one based on the energy conserva- 
tion law. Having in mind relation ([TV]), we could say that the obtained power-like 

2 Recent analysis of mean multiplicity with power-low growth in Au+Au collisions at RHIC is 
given in [ |l7| | 
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dependence which takes into account unitarity effects could be considered as a 
kind of a saturated upper bound for the mean multiplicity growth. 
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